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(54) Title: METHOD FOR MEASURING THE WAVE ABERRATIONS OF THE EYE 



Probing the eye with 
the laser beam 



Detectioa of the radiation 
scattered by retina 



Measurement of the wave front tilt in a discrete set of pupil points: 
• Discrete set of points is located alon^ concentric circles fClaim 



or in parallel to one of the axes o f a rectangular system of coordinates fClaim 4i 



< 



O 



Wave front approximation based on functions of pupil coordinates: 

Determination of oartial derivatives in anv pupil point is performed bv means of 
spline approximation fClaim li according to mathematical formula fClaim 6): 

The wave front is reconstructed usine numerical integratlort alone the radii with 
initial integration point in the center of the pupil fClaim li according to 
mathematical formula fClaim 7): 

Approximation of partial derivatives at the beeinnin^ is performed alont? each 
circle and then alone the radii fClaim 3>. or alone each of the parallel lines and 
then in orthoeonal directions fClaim 5) 



Calculation of the wave aberrations 

Based on wave ^nt reconstruction in the 
form of splines fClaim 8i 



(57) Abstract: The invention relates to medical 
instrumentation, in particular to diagnostic mea- 
suring devices, and can be applied, for exam- 
ple, for high-accuracy vision correction. This 
method uses probing the eye with a thin laser 
beam, detecting the radiation scattered by retina, 
measuring the wave front tilt in the form of first 
partial derivatives along the coordinates in a dis- 
crete set of the pupil points with known coordi- 
nates, approximating the wave front from said 
data as functions of pupil coordinates and cal- 
culating the wave aberrations of the eye. Par- 
tial derivatives are determined in any point of 
the pupil by means of spline approximation us- 
ing the values in a discrete number of points, in 
which the wave front tilts are measured. This set 
of points can be located along concentric circles 
or along one of the orthogonal axes. The wave 
front is reconstructed using numerical integra- 
tion along the radii with the initial integration 
point in the center of the pupil. The wave aber- 
rations are calculated on base of wave front data 
reconstructed in the form of splines. 
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METHOD FOR MEASURING 
THE WAVE ABERRATIONS OF THE EYE 

BACKGROUND OF THE INVENTION 

The invention relates to medical instrumentation, in particular to diagnostic measuring 
devices, and can be applied, for example, for high-accuracy vision correction. 

Methods are known for measuring the wave aberrations of the optical system of the 
human eye as a function of spatial pupil coordinates. Said methods use probing the eye with a 
thin laser beam, its backscattering by the retina, detection of the component exitmg from the 
eye by means of photosensitive devices, and measurement of the wave front tilt in different 
pupil points with known coordinates. These data are used for the wave front approximation 
and calculation of the wave aberration of the eye as wave front deformations. 

Some methods use measurement of the wave front structure at the exit of the eye 
sunultaneously in time by means of partitioning this structure into subapertures (D. R. 
Williams, et al. Rapid, automatic measurement of the eye's wave aberrations. US Patent 
6199986. Int, CI A61B 3/10, 13,03,2001). Subapertures are formed using a matrix of 
coaxial lenses and a matrix of position sensing photodetectors installed in their foci. The 
wave front is reconstructed from the measured set of tilts simultaneously in all subapertures. 

Ray tracing method for measuring the wave front and refraction aberrations 
successively in time (V. V. Molebny, et al. Device for measuring refraction aberrations of the 
eye. (Ukrainian Patent Application No. 98105286, now Patent No. 46833, Int. CI. A61B 
3/00, A61B 3/10, A61B 3/14, filed Oct. 7, 1998), is also known. According to this method, 
the eye is probed by a thin laser beam (its cross-section at the eye entrance is 0,2-0,3 mm), 
shifted in parallel to itself over the entrance aperture of the eye successively in time, the 
coordinates of the points on the retina, to which the laser beam is projected, are measured at 
each probing. Said data obtained in a discrete set of points of the eye aperture are used to 
calculate the wave front tilt in these points and then approximate the entire wave front 
surface. 

Both of the above mentioned methods use the approximation of the wave front using a 
series of Zemike polynomials. The use of said series results in some restrictions, which lead 
to non-xmiform spatial resolution when describing the refraction non-homogeneity over the 
eye aperture. Thus, in the central zone, which is the most significant for the acuity of vision. 
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spatial resolution is lower than in the peripheral zone if describing the refraction non- 
homogeneity by means of Zemike polynomials. 

To obtain a higher spatial resolution in the center, the number of polynomials must be 
increased, therefore, measurements in the larger number of aperture points are required, not 
only in the center but also at the periphery, maJdng the procedure of measurement more 
complicated: radiation exposure of the eye is increased, calculation process becomes 
significantly longer and more complicated (the number of equations increases to a few 
hvmdreds). 

Both above mentioned methods are equivalent from the point of view of wave front 
approximation using the series of Zemike polynomials. Each method can be regarded as a 
prototype. Let us take the ray tracing principle for a prototype. 

SUMMARY OF THE INVENTION 

A goal of this invention is formulated as imiformly increasing the spatial resolution of 
the measurement of the refraction non-homogeneity over the entire aperture of the eye not 
entering into a complicated embodiment of said method that could disproportionately 
increase the nvimber of probing points in the peripheral zone, increase the time of 
measurement, as well as increase the radiation dose of the patient's eye. 

The formulated goal is achieved in such a way that in the method based on the 
probing the eye with a thin laser beam, on the detection of the radiation scattered by the 
retina, on the measurement of the wave front tilt in the form of the first partial derivatives 
along the coordinates in a discrete set of pupil points with known coordinates, on the 
approximation of the wave front in the form of fimctions of pupil coordinates using said data, 
and on the calculation of wave aberrations of the eye as deformations of approximated wave 
front, the partial derivatives at any pupil point are determined by means of spline 
approximation using the values in a discrete set of points, where the wave front tilts are 
measured, and the wave front is reconstructed by means of the numerical integration along 
the radii with the initial point of integration located in the center of the pupil, where identical 
initial values of the integral are taken for all radii. 

Furthermore, a discrete set of points, in which the wave front tilt is measured, is 
located along concentric curcles with the common center, which coincides with the center of 
coordinates. 
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Fxirthermore, the approximation of partial derivatives at the beginning is performed 
along each concentric circle, and then - along the radii. 

In anotiier version, a discrete set of points, at which the wave front tilt is measured, is 
located along straight lines which are parallel to one of the axes of the rectangular system of 
coordinates. 

Approximation of the partial derivatives is performed along each of the parallel lines 
and then - in orthogonal directions. 

In addition, in still another embodiment of the proposed method, spline approximation 
along each coordinate is performed in accordance with the formula: 



(^1+1 -^i) 

where S{x) is the interpolation cubic spline along the generalized coordinate x; x/, jch-i are 
the coordinates of the pupil points i and (/ + 1) from their discrete set, at which the wave front 
tih is measured having the corresponding values S{x^), iS(x,+i); and 5'(x,.), S'ix^^i) are the 
values of the first derivative in the points x,, x^+i, which ensure the continuity of the second 
derivative S''{x) in these points. 

Furthermore, numerical integration is proposed to be performed in accordance with 
the formula: 

fr(p, <D) = FF(o, o)+ jMao) 

0 9P 

where W{p^(p) is the wave front ftmction, p - the coordinate along the radius, q> - the 
coordinate along the angle, (P,<1>) are the coordinates of the arbitrary pupil point, and 
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0 3p 



1=1 



N 



dW\ 




the interval [0,P] being divided into iV^ partial portions p J, so that / = . 



Spatial resolution in the central zone is increased in the proposed method by using the 
wave jfront approximation uniform on the entire eye without the complication of the 
measurement procedure as well as without the increase of the eye irradiation dose. 

Allocation of a discrete set of points along concentric circles with the common center, 
which coincides with the center of coordinates, helps to simplify the process of 
approximation of the partial derivatives (at the beginning, along each circle and then - along 
the radii). 

In still another version of the embodiment^ the process of approximation of partial 
derivatives is simplified due to the allocation of a discrete set of points along straight lines 
which are parallel to one of the axes of the rectangular system of coordinates and by 
performing the approximation of partial derivatives at the beginning along each of the 
parallel lines and then - in orthogonal directions. 

Uniformity of spatial resolution over the entire aperture is ensured by using the 
interpolation cubic spline S{x) in accordance with the above mentioned mathematical 
formula and under the condition of continuity of its second derivative iS"'(jc) in the 
subaperture points, where the wave front tilt is measured. 

Important peculiarity of the proposed method is also tiie procedure of numerical 
integration while reconstructing the wave front in accordance with the above mentioned 
mathematical formula along the radii from the center of the pupil in such a way that the initial 
value of integral is the same for all radii. It essentially facilitates the problem of the 
correctness of the results in the central zone, occurred if using the approximation by means of 
Zemike polynomials. 
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BIOEF DESCRIPTION OF THE DRAWINGS 

The subject matter of the invention is illustrated by drawings including 17 figures. 

Fig. 1. Course of rays in the device for measuring aberrations using Hartmaim-Shack 
sensors. 

Fig. 2. Course of rays in the ray tracing device for measuring aberration in successive 
moments of time. 

Fig. 3. Set of points for measuring the wave firont tilt located on straight lines parallel 
at least to one of the axes of coordinates. 

Fig. 4. Set of points for measuring the wave firont tilt located on concentric circles. 

Fig. 5. Measurement of the wave firont tilt using the coordinates of the laser beam 
projection on the photodetector. 

Fig. 6. Sequence of procedures while measuring the wave aberrations of the eye. 

Fig. 7. Circular trajectory of approximation. 

Fig. 8. Approximation result (broken curve) along the circular trajectory by means of 
the least squares method. 

Fig. 9. Retinal spot diagram (for the ray tracing method with the trajectory, shown in 
Fig. 7) after approximation by means of the least squares method (Zemike polynomials are 
used) - broken curve; and after approximation using bicubic splines - continuous curve. 

Fig. 10. Retinal spot diagram (for the ray tracing method with four circular 
trajectories at the entrance aperture of the eye) after approximation by means of the least 
squares method (using Zemike polynomials). 

Fig. 11. Retinal spot diagram (for ray tracing method with four circular trajectories at 
the entrance apertm-e of the eye) after approximation using bicubic splines (for the same 
initial data as shown in Fig. 10). 

Fig. 12. Sequence of procedures while measuring the wave aberrations of the eye 
using the proposed method. 

Fig. 13. Functional schematic diagram of a device for measurement of the wave 
aberrations of the eye using the proposed method. 
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Fig. 14. Schematic design of the optical block of device for measuring the wave 
aberrations of the eye using the proposed method. 

Fig, 15. General view of a device for measuring the wave aberrations of the eye using 
the proposed method. 

Fig. 16. Map of refraction errors, reconstructed by means of the least squares method 
using Zemike polynomials. 

Fig. 17, Map of refraction errors, reconstructed when using the wave front 
approximation with the bicubic splines. 

Let us discuss in detail implementation of the proposed method for measxiring the 
wave aberrations of the human eye. 

Fig. 1 explains the concq)t of measurement of the wave aberrations of the eye using 
Hartmann-Shack sensor. The optical system of the eye is represented in this figure as a thin 
lens and the retina is placed in the focus of this lens. A thin parallel laser beam is directed 
into the eye. Preferably the beam enters the central zone of the pupil so that its axis coincides 
with Ihe visual axis of the eye (for simplicity, physiological peculiarities of the vision are not 
taken into account, according to which, visual and optical axes are to be differentiated - we 
shall use these terms as synonyms). The optical system of the eye focuses the laser beam into 
the point whose size is determined in an ideal case only by the refraction. 

This illiuninated point acts as a secondary point irradiator whose light is scattered in a 
certain spatial angle and a part of it is propagatuxg through the eye piq)il. A matrix (raster) of 
coaxial lenses with diameter approximately 0,3-1,0 mm (sometimes even less) is installed on 
the path of radiation. These lenses have identical focal lengths. System of position sensing 
photodetectors is installed in the foci of the lenses. It is common to use a TV matrix of 
charge-coupled devices (CCD). 

Coordinates of each illuminated point are measured in foci of microlenses. Deviation 
of the illuminated point from the optical axis of a corresponding lens is the measure of the 
deviation of the wave front in the zone of this lens (in the subapertute of the wave front), 
being (in linear approximation) the first partial derivative along the coordinates for this wave 
front subaperture. Based on these data (deviation and coordinates), the wave front is 
approximated as fimctions of pupil coordinates and then aberrations of the wave front are 
calculated as its deformations. 
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Ray tracing method differs in several details, but the sequence of procedxires is the 
same as with Hartmann-Shack sensors. The main difference is successive (in time) 
measurement of the wave firont tilts in its subapertures. 

Fig. 2 shows the course of rays in the ray tracing procedure in successive moments of 
time. In this figure, flie optical system of the eye is shown as a thin lens, the retina being 
located in its focus. In the first moment, a thin laser beam is directed into the eye at the pupil 
point with coordinates (xi, y{). In the next moment, said beam is directed to a point with the 
coordmates (X2,y2\ The procedure continues in this manner up to the last point of the eye 
aperture with the coordinates Qcn, y^). Each tune, the laser beam is projected on the retina as a 
point (a point is an abstraction, actually, it is a more sophisticated figure, though roughly 
speaking it may be considered as a point). 

Like in the above case, this point is a secondary point irradiator, and it is possible to 
detect it by means of a photodetector with an objective lens. Also, it is a position sensing 
photodetector generating at its output tibie coordinates of the image of the retinal point, into 
which the laser beam is projected after having entered the eye pupil with the known 
coordinates (jc,, yi). 

Deviation of the laser beam projection dc on the retina is unambiguously related to the 
measured deviation in the plane of photodetector. This value is the measure of beam 
deviation or adequately the measure of wave firont deviation at the point of eye aperture with 
known coordinates. Thus, like in the first method, as a result of performed measurements, a 
wave firont tilt is obtained in the points of eye aperture with the known coordinates. 

Based on these data (like in the first method, this is the tilt and the coordinates), the 
wave front is approximated as fimctions of pupil coordinates, and then aberrations of the 
wave front are calculated as its deformations. 

Figs. 3 and 4 depict the disposition of the set of points on the eye aperture, where the 
wave front tilt is measured, the points of the set in Fig. 3 being located along straight lines in 
parallel to the axes of the Cartesian system of coordinates. Their coordinates (x/, y{) are also 
represented in the rectangular system of coordinates. In Fig. 4, these points are located along 
the concentric circles. It is convenient to describe their coordinates in the polar system of 
coordinates (/?/, <pi). 

As a rule, the rectangular system of coordinates is used in the method based on 
Hartmann-Shack sensor (Fig. 3), and the polar system of coordinate, more appropriate for a 
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natural round eye aperture, is used in the ray tracing method (Fig. 4). The use of one or 
another system of coordinates in the proposed method affects only the procedure of 
calculation of the approximation functions. 

Fig. 5 shows the relation between the measured value of the shift of the laser beam 
projection in a subaperture (Hartmaim-Shack sensor method) or in the entire aperture (ray 
tracing method). The wave front tilt can be described by an angle a or by a imit vector n 
normal to the wave front surface. The angle a is determined from the relation for a triangle 
ODO': 

O'D 

where 0'D=8^ is the measured value of the shift of the laser beam projection in the 

photodetector plane, and the focal length of the objective lens. One gets from this 

expression: 

The sequence of the procediares in the course of measurement of wave aberrations is 
illustrated in Fig. 6. This sequence comprises the following steps: 

(a) probing the eye with the laser beam; 

(b) detection of the radiation scattered by the retina; 

(c) measxirement of the wave front tilt in a discrete set of pupil points; 

(d) approximation of the wave front with the ftmctions of pupil coordinates; 

(e) calculation of wave aberrations. 

In the prototype method in order to approximate the wave front, series of Zemike 
polynomials are used. The solution of the problem is assisted by the least squares method. 
This procedure is schematically illustrated in Figs. 7 and 8. A circular trajectory of scanning, 
which means the trajectory of constant radius p, is shown in Fig. 7. This trajectory includes 
fixed points with known coordinates (p, tpi), (p, <p2\ ip, 9d> ip. <Pn\ in which the value 
5^ is measured. 

Diagram of values d^, for points from <p\ to (p^ is shown in Fig. 8. A module of this 
value is shown in the diagram, its orthogonal components being analyzed in the analytical 
record. A continuous ftmction (series of polynomials) is illustrated in Fig. 8 by dotted hne 
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which, by means of the least squares method, approximates the trajectory of beam projection 
on the photodetector (for the ray tracing method). As one can see, the approximation function 
takes the values differing jfrom the measured values. This is one of the sources of the error of 
measurement. Another source results from the use of 2femike polynomials leading to a non- 
uniform spatial resolution when describmg refraction non-homogeneity in the eye aperture, 
whereas in the central zone the spatial resolution of measurements is lower than in the 
peripheral zone, that reduces the valuables of the obtained information on eye aberrations 
since this central zone is the most significant for the acuity of vision. 

A higher spatial resolution can be reached by increasing the number of measuring 
points, allowing a larger number of approximating terms of the series. This will usually lead 
to more complicated measurements and will increase the dose of retina irradiation. The 
calculation process also becomes more complicated and longer due to necessary increase of 
the number of equations. 

Spline approximation is used in the proposed method, according to which, the 
trajectory of beam projection on the retina is not approximated by means of the least squares 
method (broken trajectory Fig. 9), but passes through the measured values of beam deviation 
for probmg points with coordinates (p, q>i\ {p, (pi), (p, ^,.), (9, ^;^), 

Figs. 10 and 11 are examples of trajectory of beam projection on the retina (retinal 
spot diagram) approximated using the same measurement data: by means of the least squares 
method in Fig. 10, and by means of splines in Fig. 11. Flatness of the approximated 
trajectories in Fig. 10 is obviously more distinct than in Fig. 11, that certifies the loss of 
considerable part of information on variations of the refiraction non-homogeneity containing 
the higher spatial frequencies. 

The proposed method for measurement of wave aberration of the eye differs from its 
prototype due to inclusion of some distmctive features (as illustrated in Fig. 12) into the 
procedures shown in Fig, 6, which are essential for the proposed method. Let us discuss 
aforementioned distuictions in the sequence of their use referencing to the corresponding 
claims. 

In the procedure of measurement of the wave front tilt in a discrete set of the pupil 
points, this set of points is formed either along concentric circles (as provided in Claim 2) or 
in parallel to one of the axes of the rectangular coordinate system (Claim 4). 

9 
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The wave front approximation using functions of pupil coordinates includes the 
following distinctive features: 

(a) Partial derivatives in any point of the pupil are determined using spline 
approximation (Claim 1). Accordmg to the Claim 6, this approximation is performed 
according to the relation, in which the interpolation cubic spline S{x) along the generalized 
coordinate x is derived from the data on the coordinates x„ x^-i of the pupil points i and (/+1) 
fix)m their discrete set, in which the wave front is measured, and from the data on the values 
of the first derivative in the points x,, x^. The condition of continuity of its 

second derivative S''(x) in these points must be satisfied. 

(h) Wave front is reconstructed by means of the numerical integration along the radii 
with the initial integration point in the pupil center (Claim 1). The Claim 7 specifies this 
procedure by means of mathematical relation for the wave front W(P,0) in any point of the 

pupa usmg known values of derivatives ^^^'^) in the polar system of coordinates {p,<p). 

The interval of integration [0,P] is divided here into impartial portions |>,_i, p J so that 

/ = 1, , and the intergal f^^'^^ dp is replaced by a sum 

0 

(c) Approximation of partial derivatives, at the beginning, is performed along each 
circle and then - along the radii (Claim 3). In another embodiment of this method, the 
approximation is performed along each of the paraUel straight lines and then - in orthogonal 
directions (Claim 5). 

hi the process of calculation of wave aberrations, an essential distinctive feature is 
also included: wave aberrations, for example, in the form of a series of Zemike polynomials 
are now calculated not from a limited set of data including only the discrete points, in which 
the measurement was taken, but the continuous fimctions. Thus, for calculation of any 
number of the terms of a Zemike series, there is no need in increasing the number of 
measuring points, since the wave front surface is described by splines complying with the 
continuity conditions. 
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A functional schematic diagram of the device for measurement of wave aberrations of 
the eye based on this method is shown in Fig. 13. It comprises: a probing channel 1, a 
detecting channel 2, a channel 3 of information processing and displaying, a channel 4 of 
orientation of the patient's eye, and a polarization beam splitter 5 installed at the exit of the 
probing channel 1, and a first beam splitter 6 installed in firont of said detecting channel 2 and 
said channel 4 of orientation of patient's eye. An eye imder investigation 7 comprises a 
cornea 8, a crystalline lens 9, and a retina 10. 

Said probing channel 1 consists of successively installed laser 1 1 and a imit of beam 
positioning 12. Said detecting channel 2 includes successively installed an optical unit 13 and 
aphctodetector 14. 

Said chaimel 3 of information processing and displaying is made of a unit 15 of spline 
approximation, a unit 16 of wave firont reconstruction, a unit 17 of calculation of wave firont 
aberrations, a unit 18 of information display, and a unit 19 of coordinate data. 

The input of said imit 15 of spline approximation is connected to the output of said 
photodetector 14. The units 15, 16, 17 and 18 are connected in such a way that the input of 
each next unit is connected to the output of the previous one. The output of said unit 19 of 
coordinate data is connected to said unit 12 of beam positioning and the input of said unit 15 
of spline approximation. 

Said chaimel 4 of orientation of patient's eye consists of a second beam splitter 20, a 
TV camera 21, a light source 22 and an accommodation transparency 23, here the elements 
22, 23 are installed successively, and the radiation of the light source is directed into the 
patient's eye so that it is reflected twice: by said second beam splitter 20 and by said first 
beam splitter 6, and passes through said polarization beam splitter 5. 

Functioning of the device. Before measurement, said patient's eye 7 is to be oriented 
so that its line of sight coincides with the optical axis of the device. For this purpose, the 
patient has to direct his/her sight onto the light source 22 which is made of a Ught source 
itself which can be, for example, a light emitting diode (LED), and of an optical forming 
system. When using this system, the accommodation transparency 23 can be conjugated with 
the retina of the eye in a wide range of distances (for the sake of simplicity, these elements 
are not shown in Fig. 13). Said TV system 21 ensures the visualization of the orientation 
process. This process is important for obtaining correct results of measurements but not 
essential for showing distinctive features of the proposed method. 
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In the moment when the patient's eye is correctly oriented and accommodated to the 
target distance, the laser 11 is switched on. Laser radiation is stepwise shifted all over flie 
pupil as shown in Figs. 2 and 4. After passing through the refracting surfaces and eye 
structures 7, this radiation hits the retina. Radiation of the laser 1 1 is linearly polarized. As a 
result of scattering in the retinal structures, this radiation is depolarized and only the cross 
component passes through the polarization beam splitter 5 on the way back to the detecting 
channel 2. 

A part of laser radiation reflected by the comeal surface 8 does not pass throu^ the 
polarization beam splitter 5 (it is well known that the linear polarization is not changed after 
such reflection). Therefore, only the cross component of radiation scattered by the retina 10 
reaches the photodetector. 

The plane of the retina 10 and the plane of the photodetector 14 are conjugated by 
means of said optical unit 13. The measured values of deviations (Figs. 2 and 5) are sent to 
said unit 15 of spline approximation. Actually, the unit 19 of coordinate data is the driver of 
said unit 12 of beam positioning which controls the spatial position of the laser beam that 
enters the eye in parallel to the optical axis of the device. Besides, the data on the coordinates 
of the pupil points in which the wave front tilt (partial derivatives) is measured in this 
moment, are transmitted from said unit 19 of coordinate data to said unit 15 of spline 
approximation. The functions of pupil coordinates, approximating the values of the partial 
derivatives (unambiguously defibtied from the deviations in any point of the pupil, are 
reconstructed in said unit 15. 

Said unit 16 of wave front reconstruction performs its function using numerical 
integration along the radii with the initial integration point in the center of the pupil, where 
the identical initial values of the integral are defined for all radii. 

Said unit 17 of calculation of wave front aberrations analyses the deviations of the 
wave front from an ideal form, and the map (two-dimensional distribution) of said deviations 
is visualized by said unit 18 of information display. 

The structure of the above described device and interaction of its main components 
are based on the ray tracing method. The stmcture of a device according to the method which 
uses the partitioning of the radiation, exiting from the eye, into subapertures, would be in 
general the same, only said optical unit 13 should perform functions of Hartmann-Shack 
sensor. Said unit 12 of beam positioning would perform the functions of beam shaping, and 
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the information on the coordinates of the subapertures for said unit 19 of coordinate data 
should be formed outside of said unit in the detecting channel 2. 

Any modification declared in the dependent claims can be implemented using the 
described device. 

For example, said unit 12 of beam positioning performs the function of allocation of a 
discrete set of points, at which the wave front tilt is measured, along the concentric circles 
(see Fig. 4) with the common center, which coincides with the center of coordinates (Claim 

2) . SpUne approximation of partial derivatives at the beginning is performed along the 
concentric circles using said unit 15, and then - along the radii (Claim 3). 

When implementing the proposed method with the Hartmann-Shack sensor, a discrete 
set of points to measure the tilt of the wave front is located along the straight lines (see Fig. 

3) , being parallel to one of the axes of the rectangular system of coordinates. Such 
positioning, as a rule, corresponds to the structure of the matrix raster of said optical xmit 13. 
This structure of a set of the points ensures said unit 15 of spline approximation to perform 
the approximation of the partial derivatives along each of the parallel lines - at the begiiming, 
and after that - in orthogonal directions (Claim 5). 

It is necessary to stress, that the straight structure of a set of points can also be 
realized in the ray tracing method but the circular structure is more rational for it. 

Spline approximation along each coordinate is performed using said unit 15 in 
accordance with the formula: 

, (xj+l-xfix-Xj) r ^ (izfiZi^illJliO eYv ^ 
-XiP {Xi+i -XiY 

where S{x) is the interpolation cubic spline along the generalized coordinate x; x/, Xi^y are 
the coordinates of the pupil points / and (i + 1) from their discrete set, in which the wave 
front tilt is measured having corresponding values 5(x,), S{xm); and 5"(x,), are 
the values of the first derivative in the points x„ xh-i, which ensure continuity of the second 
derivative S"{x) in these points (Claim 6). 
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Wave front reconstruction using said unit 16 can be performed by means of numerical 
integration in accordance with the formula: 

^r(p,<i») = wio,o)+ /Mao) jp. 

0 

where W(p,((>) is the wave front function, p - the coordinate along the radius, (p - the 
coordinate along the angle, (P.O) are the coordinates of an arbitrary pupil point, and 

0 '^P f=l 

the interval [0,P] being divided into impartial portions p,] so that / = ijif (Claim 7). 

Finally, said unit 17 of calculation of wave aberrations can calculate the wave 
aberrations not from the data for limited number of points of flie pupil, at which the wave 
front tilts are measured, but from the wave front data reconstructed in the form of splines. 
That is why there is no need in increasing the number of probing points and consequently in 
increasing the radiation exposure of the patient's eye, that ensures indisputable advantage of 
the proposed method together with the increase of the spatial resolution of the measurement 
of the wave aberrations. 

Fig. 14 illustrates an example of a construction of said channels 1, 2 and 4 of the 
device, hnplementing the proposed method. All components of the above mentioned channels 
are installed on a rigid plate fixed to a movable base 5 of a slit-lamp. Patient has to take place 
to the right side of the construction having a chin rest CR and a forehead rest FR, to which the 
patient presses his chin and his forehead. 

The construction has three degrees of freedom: up-down, left-right, and forward- 
backward. Movements and turnings of the patient's head help to reach the comcidence of the 
optical axis of the device with the line of sight of the eye under investigation, and then the 
measurements are made. A complete cycle of probing the eye in a complete set of the probing 
points lasts few tens of miUiseconds. 

A general view of the device is shown in Fig. 15. Besides the slit-lamp base (Fig. 14), 
the device also comprises a PC-based channel of information processing and displaying. Not 
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only the wave front map but also additional information can be demonstrated on the PC 
display, this information can be overlaid/combined with the output information from the TV 
system 21. 

The information on the wave front aberrations can be converted into more easily 
perceptible information on the refraction errors. Examples of the information presented in the 
form of the two-dimensional color-coded maps are shown in Figs. 16 and 17. In Fig. 16, the 
map is reconstructed as a result of wave front approximation using Zemike polynomials by 
means of the least squares method, and in Fig. 17, the map is reconstructed from the same 
primary measurement data but is the result of spline approximation. Comparing said maps 
makes it obvious that in the first case, the very important information on keratoconus is 
missing, but is clearly seen in Fig. 17 in the left lower quadrant. This example proves evident 
advantages of the use of the proposed method for measurement of wave aberrations of the 
eye. 
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CLAIMS 

1 . A method for measuring the wave aberrations of the eye, based on probing the 
eye with a narrow laser beam, detection of the radiation scattered by the retina, measurement 
of the wave front tilt in the form of the first partial derivatives along the coordinates in a 
discrete set of pupil points with known coordinates, approximation of the wave front in the 
form of functions of pupil coordinates using said data, and calculation of wave aberrations of 
the eye as deformations of the approximated wave front, characterized 'm that the 
partial derivatives at any pupil point are determined by means of spline approximation using 
the values in a discrete set of points, where the wave front tilts are measured, and the wave 
front is reconstmcted by means of the numerical integration. 

2. A method according to claim 1, wherein the wave front reconstructed by 
means of the numerical integration is along the radii with the initial point of integration 
located in the center of the pupil, where identical initial values of the integral are taken for all 
radii. 

3. A method according to claim 1, characterized 'm that the discrete set 
of points, in which the wave front tilt is measured, is located along concentric circles with the 
common center, which coincides with the center of coordinates. 

4. A method as in claim 3, characterized in that the approximation of 
partial derivatives at the beginning is performed along each concentric circle and then along 
radii. 

5. A method according to claim 1, characterized in that the discrete set 
of points, in which the wave front tilt is measured, is located along straight lines, which are 
parallel to one of the axes of the rectangular system of coordinates. 

6. A method as in claim 5, character izedin that the approximation of the 
partial derivatives is performed along each of the parallel straight lines, and then in the 
orthogonal directions. 

7. A method as in claims 1, or 4, or 6, c A a racier iz^rfin that the spline 
approximation along each coordinate is performed in accordance with the formula: 
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(^1+1 -^/T (^1+1-^/) 

where ^(x) is the interpolation cubic spline along the generalized coordinate x; x/, xm are the 
coordinates of the pupil points i and + 1) from their discrete set, at which the wave front tilt 
is measured having corresponding values S{xf^^); and iS'(jc,), 5'(x,.+i) are the values 

of the first derivative in the points x/, x^-u which ensure continuity of the second derivative 
S''{x) in these points. 

8. A method according to claim 2, characterized in that the numerical 
integration is performed in accordance with the formula: 

»r(P.«i>)=Fr(o.o)+jM^^p, 

where W{p,i^) is the wave front function, p - the coordinate along the radius, 9 - the 
coordinate along the angle, (P, O) are the coordinates of arbitrary pupil point, and 

^ fdWjp,^) ^ ^ L 2 ' )f . . 

J ^Z-^<^P « S ^"—J- -{Pi - 

0 °P f=i °P 

the interval [0,P] being divided into iV^ partial portions pj, so that i = . 

9. A method according to claims 1, 2, 3, 4, 5, 6 or 8, characterized in 
that the calculation of wave aberrations is performed based on the wave front data 
reconstmcted in the form of splines. 

10. A method according to claim 7, characterized in that the calculation 
of wave aberrations is performed based on tiie wave front data reconstmcted in the form of 
splines. 
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11. A method according to claim 1, characterized m that the reconstructed 
wave front is used for the purposes of vision correction. 

12. A method for measuring the wave aberrations of the eye, based on probing the 
eye at a discrete set of pupil points with a narrow laser beam, detection of the radiation 
scattered by the retina, measurement of the wave front tilt values in the form of first partial 
derivatives along known coordinates of the discrete set of pupil points in a coordinate system, 
approximation of the wave front in the form of functions of pupil coordinates using data from 
the detection of radiation scattered by the retina, and calculation of wave aberrations of the 
eye as deformations of an approximated wave front, wherein the first partial derivatives at 
any pupil point are determined by means of spline approximation using the data for the 
discrete set of pupil points, where the wave front tilt values are measured, and the wave front 
is teconstructed by means of numerical integration. 

13. A method according to claim 12, wherein the wave front reconstructed by 
means of numerical integration is in a polar coordinate system and is along the radii to the 
pupil points with the initial point of integration located in the center of the pupil, where 
identical initial values of the integral are taken for edl radii. 

14. A method according to claim 12, wherein the discrete set of pupil points, in 
which the wave front tilt is measured, is located in a polar coordinate system along concentric 
circles with a conmion center, which coincides with a center of the radial coordinates. 

15. A method as in claim 14, wherein the approximation of partial derivatives at 
the begmning is performed along each concentric circle and then along radii. 

16. A method according to claim 12, wherein the discrete set of pupil points, in 
which the wave front tilt is measured, is located in a rectangular system of coordinates along 
straight lines, which are parallel to one of the axes of the rectangular system of coordinates. 

17. A method as in claim 16, wherein an approximation of the partial derivatives 
is performed along each of the parallel straight lines, and then in orthogonal directions. 

18 A method as in claims 12, 15 or 17 wherein the spline approximation along 
each coordinate is performed in accordance with the formiila: 
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where iS'(x) is the interpolation cubic spline along the generalized coordinate jc; x^, xm are the 
coordinates of the pupil points / and (/ + 1) from their discrete set, at which the wave front tilt 
is measured having corresponding values S{xf^^); and 5"(jc,.), are the values 

of the first derivative in the points x,, xm, which ensure continuity of the second derivative 
iS''(x) in these points. 

19. A method according to claim 13 wherein numerical integration is performed in 
accordance with the formula: 

0 

where W(p,(p) is the wave front function, p - the coordinate along the radius, (p - the 
coordinate along the angle, (P,<I>) are the coordinates of arbitrary pupil point, and 



0 M °P 



the interval [0,P] being divided into i\r partial portions p,], so that i^X^N . 

20. A method according to claims 12, 13, 14, 15, 16, 17 or 19, wherein the 
calculation of wave aberrations is performed based on the wave front data reconstructed in 
the form of splines. 

21. A method according to claim 18, wherein the calculation of wave aberrations 
is performed based on the wave front data reconstructed in the form of splmes. 

22. A method according to claim 12, wherein the reconstmcted wave front is used 
for the purposes of vision correction. 
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Probing the eye with 
the laser beam 



Detection of the radiation 
scattered by retina 



Measurement of the wave front tilt in a discrete set of pupil points 
• Discrete set of p oints is located alons concentric circles 



(Claim 2) 



or in varallel to one of the axes o f a rectangular system nf coordinates (Claim 4) 
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spline approximation (Claim 1) according to mathematical formula (Claim 6): 
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